Network connectivity through small openings 
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Network connectivity is usually addressed for convex domains where a direct line of sight exists 
between any two transmitting/receiving nodes. Here, we develop a general theory for the network 
connectivity properties across a small opening, rendering the domain essentially non-convex. Our 
analytic approach can go only so far as we encounter what is referred to in statistical physics 
as quenched disorder making the problem non-trivial. We confirm our theory through computer 
simulations, obtain leading order approximations and discuss possible extensions and applications. 



I. INTRODUCTION 

Wireless ad hoc and sensor networks [1] have attracted 
much attention in recent years as they do not rely on a 
pre-existing system infrastructure, such as central routers 
or access points. Instead communication devices, here 
referred to as nodes, can relay messages to one another 
in a multi-hop fashion, thus achieving good coverage and 
connectivity over large areas. Modern applications of 
wireless ad hoc and sensor mesh networks can be found 
inter alia, in swarm robotics [2] deployed in hazardous 
locations such as factories, mines and disaster areas [3], 
laptops, power distribution ("smart grid") technologies 
[4], vehicles for road safety [5], and underwater sensor 
networks [6]. 

Achieving smooth functionality of a network in the 
form of good connectivity, whilst simultaneously mini- 
mizing economic, computational and energy costs is a 
difficult task, addressed from a physical layer's perspec- 
tive by the theory of random graphs [7] and complex net- 
works [8]. Here, networks are modelled by a collection of 
nodes randomly distributed in a region of space, pair- 
wise connected with a relative position dependent prob- 
ability. The network's connectivity properties can then 
be characterized via a plethora of methods and measures 
such as for example various clustering statistics, network 
modularity measures, node importance and correlations 
between degrees of neighbouring nodes. Of particular in- 
terest in dense networks is the optimal number of nodes 
necessary to maintain full connectivity with certain prob- 
ability [9]. 

A theory for predicting the probability Pfc with which 
a random network is fully connected in confined convex 
spaces was recently developed in [10] and [11]. There, it 
was shown that Pfc is strongly influenced by the details 
of the domain boundary and in fact is asymptotically 
dominated by the lack of connectivity near the sharpest 
interior corner. The main result was a general analytic 
formula facilitating the design and analysis of reliable 
wireless mesh networks [12]. 

In this paper, we extend this theory to non-convex do- 
mains and in particular to domains involving a small 
keyhole type opening as shown in Fig. 1. The notion 
of the keyhole problem in wireless communications was 
introduced in [13], and has since been extensively studied 
from the antenna propagation perspective and the rele- 




FIG. 1: Clustering of a random network formed in a square 
domain of side L = 10 at density p = 2. A hard vertical wall 
with a keyhole opening of width w = 1 splits the domain into 
two equal parts thus hindering full connectivity. 



vant channel distributions [14]. Instead, here we focus 
on the network aspect of the keyhole problem and ex- 
plain why this system exhibits semi-quenched disorder, 
a feature common in spin glass systems (e.g. disordered 
magnets with frustrated interactions) [15]. As our main 
contribution, we derive a general formula for the prob- 
ability of two sub-networks to connect through a small 
opening. We confirm our results using computer simu- 
lations and furthermore obtain leading order asymptotic 
expansions in closed form. To the best of our knowledge, 
together with the recent works of [16], our results consti- 
tute the first analytic approach to network connectivity 
in non-convex domains. 

The paper is structured as follows: In Sec. II we de- 
scribe the keyhole problem set-up and define the relevant 
parameters and observables. In Sec. Ill we recall how the 
probability of full connectivity is calculated in convex do- 
mains and in Sec. IV we reduce the problem of full net- 
work connectivity to that of the existence probability of a 
bridging link between sub-domains. In Sec. V we obtain 
leading order analytic approximations to this probability 



and in Sec. VI discuss why the system exhibits quenched 
disorder and thus derive our main result, a general for- 
mula for the probability of two sub-networks to connect 
through a small opening. In Sec. VII we confirm our re- 
sults through computer simulations and discuss asymp- 
totic limitations and difficulties due to the many length 
scales associated with the problem at hand. Finally, some 
concluding remarks are given in Sec. VIII. 



II. DESCRIPTION OF THE PROBLEM 

In its simplest form, the problem set-up involves a do- 
main V C M^ of area V which is split into two convex 
sub- domains A,B C V, containing Nj\ and Njg nodes 
respectively. The two domains are separated by a hard 
straight wall containing a small keyhole opening of width 
w (see Fig. 1). For reasons to be discussed later, we as- 
sume that the keyhole is located sufficiently away from 
the remaining borders of A and B. The node locations 
a,j e A and bj e B ior i e [1,AU] and j € [1,Nb] are 
chosen at random from a uniform distribution with den- 
sity PA.B ~ Nj^^b/Va.b where V/i,e is the corresponding 
volume of the sub-domains. 

After deployment of the nodes, communication links 
between pairs of nodes are established with probability 
H{r), where r is the relative distance between the pair. 
We maintain physical relevance by adopting a specific 
pair connectedness function H{r) derived from wireless 
communication theory [17] and applicable to ad hoc and 
sensor networks - the Rayleigh fading model [22]. Here, 
the information outage probability Pout due to a single- 
input single-output (SISO) [23] system is given by Pout = 
Pr (log2(l -I- SNR X |/i|^) < a;), where x is the minimum 
outage rate threshold, SNR oc r~'^ is the signal to noise 
ratio and 77 is the path loss exponent [24]. Finally, h 
is the channel transfer coefficient and \h\'^ is modelled 
as an exponentially distributed random variable. Hence, 
the connection probability between two nodes a distance 
r apart can be expressed as 



H{r) = 1 - Po 



= e-^^\ 



(1) 



where /3 sets the characteristic connection length tq = 
(3~^/^. It is worth noting that in the limit of 77 — > 00, 
the connection between nodes is no longer probabilistic 
and converges to the well studied case in geometric graph 
theory [7] , the unit disk model with an on/off connection 
range at the limiting tq. 

Figure 1 shows a typical realization of a random net- 
work in a square domain of side L = 10 with w = 1 at 
a density of p^ = pg = 2 using /? = 1 and 77 = 00. We 
observe that two large connected components (clusters) 
are formed, one in each sub-domain, with smaller ones 
forming near some of the corners. In this paper, we will 
investigate the probability of achieving a fully connected 
network denoted by PL for a given pair of densities pj\ 
and pb- 



III. FULL CONNECTIVITY IN SUB-DOMAINS 

In [11], it was shown that the probability of a spa- 
tially confined random network being fully connected Pfc 
is dominated by local boundary effects due to "hard to 
connect to" regions. Thus, contrary to the expected uni- 
versal features of large networks, full connectivity is gov- 
erned by local and not global features of the network 
domain. Through the use of a cluster expansion, often 
used in statistical physics, it was shown that for arbitrary 
convex domain V C M"* 
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where ri2 = |ri — r2| is the relative distance between 
nodes with position coordinates ri,r2 € V. For fast de- 
caying connectivity functions H{r) (as in Eq. (1)), the 
outer integral in (2) can be approximated by a sum over 
separate boundary objects of different co-dimension i — 
0, 1, ... d (with i ~ corresponding to the bulk/volume 
term) 
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where d £ N is the space dimension, Gj. is a geometrical 
factor for each object j of co-dimension i while Vj. is the 
corresponding d— i dimensional volume of the boundary 
object with internal solid angle cjj. G (0, 27r''/^/r(d/2)). 
Equation (3), suggests that to a good approximation 
network connectivity is the sum of its parts. For example, 
using 7/ = 2 in (1) we have that a random network of N 
nodes confined in a convex n-gon of area Vi, perimeter 
V2 and interior angles fl = {u!i,u!2 . . . a;„}, will be fully 
connected with probability 



Pfc^l- Vipe' 
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where p — N/Yi. Therefore, for the keyhole set-up 
of Sec. II, it follows that full connectivity in each sub- 
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and PL is given by (4) . Notice that at high 



densities p, the term with the sharpest corner min(i7) will 
dominate Pfc- 



IV. FULL CONNECTIVITY 

If both sub-domains are fully connected, and at least 
one "bridging" link exists through the keyhole opening, 
then the combined network is also fully connected. Thus 
we may write 



p(V) ^ p(A)p{B)j. 
fc fc fc ' 



(5) 



where X is the probability of at least one bridging link 
between the two sub-domains. As with the simplification 




FIG. 2: Close-up of the keyhole opening showing a node a^ in 
sub-domain A and its associated connection region B{ai) C 
B. The dashed circle simply illustrates that the keyhole is 
sufficiently isolated from any other boundary elements of the 
domain V. 



of (2) into (3), the assumption of statistical independence 
between the three events in (5) is justified due to the fast 
decaying H{r). For w sufficiently small, no node in A is 
likely to influence the connectivity properties of nodes in 
B, and vice versa. Writing PL as in (5) is particularly 
helpful when performing computer simulations as calcu- 
lating X is typically much faster than calculating pL 
(see Sec. VII). 

Noting that the probability of at least one bridging link 
is just the complement of no bridging links, we have that 

^=i-((nn(i-x«^-^'j-))6)-4' (6) 

where Hij = H{\ai ~ hj\) and Xij equals 1 if a "line 
of sight" exists between nodes a^ and bj, and is zero 
otherwise. The angled brackets in (6) represent a spatial 
average over all possible node configurations within each 
sub-domain 

('^)C = 77]V7 / 0(Ci,C2,...,CA,^)dCidC2...dCA,(,, 

(7) 

with C — A,B, and c = a, b, respectively. Note that 
the order in which the average in (6) is performed is not 
important. 



V. LEADING ORDER APPROXIMATION 

The crudest approximation of X can be obtained by 
assuming that all spatial integrals in (6) separate such 
that 



= 1 - CXp ( -PAPB 
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where we have assumed that A'^, N^ ^ 1 and have used 



the fact that [1 — x) 



N 



-Nx 



in the limit of N 



oo. We shall be using this approximation throughout 
our analysis. 



In Eq. (8) we can identify {{xijHij)B)A as the prob- 
ability of a randomly selected pair of nodes to form a 
bridging link. We now calculate the integrals in (8) for 
a soft connectivity function H{r) with 77 = 2. In doing 
so H(r) is a Gaussian function thus rendering the math- 
ematics tractable. For every node an — {xi,yi), there is 
a line of sight connection region ^(a^) C B (see Fig. 2). 
Hence we have that jj^XijHijdhj = ^^.Jlijdhj. In 
Cartesian coordinates, we get 



B{.x,,yi) 



iJydbj = 



00 /'■U.-i- 



+ erf(v/^(y, -u_)))da; 



e-^(^'-^^)'(erf(y^(u+-y,)) 



(9) 



where u± = {yi =F w/2)xj/xi ± w/2 are as seen in Fig. 2. 
Setting the upper integration limit of dxj to infinity in 
an otherwise finite domain V, is justified by the exponen- 
tially decreasing integrand function. 

The change of variables {xi,yi,Xj) — \/y/]3{xi,yi,Xj) 
simplifies the expression in (9) whilst a further change 
of variable Zj = —Xj/xi eliminates the Xi dependence in 
the argument of the error functions. We may now change 
the order of integration and perform the Xi integral over 
the negative real line first to obtain 

erf [(1 + z,){m + ^)] - erf [(1 + %)(y, - ^)] 



2(1 + z,Y 



(10) 



which is an even function in iji. Integration in yi G 
(—(X), 00) gives \/ttw / {2{1 + Zj)'^ (3'^/'^) and finally integra- 
tion in Zj e (0,00) leads to the immense simplification 
of 



X„i?ydbjda, = -— -, 

J AJB ^P ' 



and hence that 



X « 1 - exp -pAPB 



■ 2/33/2 



(11) 



(12) 



Indeed, Eq. (12) offers a decent approximation of X at 
low densities and for w ^ 1. This is confirmed in Fig. 4 
where we compare against computer simulations of X. 
We postpone discussing the accuracy of (12) for Sec. VII. 
For a hard connectivity function H{r) with rj = 00 and 
connectivity range rg, the sector angle of B{ai) is given 
in polar coordinates by (f>i — ~wcos9i/ri + 0{w'^) (see 
Fig. 2). Hence we have that 



/ X^jH^JdhJ w -^ {nrl 
Jb 27r 



(13) 



where we have approximated B{a^) by an annular sector 
with radii r^ < tq. Integrating over A we have that 



3^/2 pro J,. 2r37;i 

/ r.'^{rl-r^)drM.= ^-^. 

7r/2 Jo ^ 'i 



(14) 
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FIG. 3: Overlapping sectors of radius ro = 1 due to nodes in A. The union of these sectors intersected with A defines 
the cumulative connection region {_]^B{ai). Note that some nodes are too far away from the keyhole opening and so do not 
contribute to U^Z?(ai). For larger openings (left) the cumulative connection region is dominated by the nearest node to the 
keyhole, while for smaller openings (middle) the cumulative connection region defines a complicated region. For very small 
openings (right), there is very little overlap between sectors. 



Notice that substituting ro = /3^^ into (14) gives the 
correct scaling as in (11). The pre-factor however 2/3 < 
■y/7r/2 w 0.886 suggests that connectivity through small 
openings is better for soft connectivity functions H[r). 



VI. SEMI-QUENCHED DISORDER 

Equation (8) suffers from the unphysical assumption 
of uncorrelated node positions. That is, if one node in B 
is situated very close to the keyhole opening, then almost 
all nodes in A will have an increased probability of con- 
necting to it and vice versa. Thus, Eq. (8) fails to capture 
this feature of the network as it assumes that the node 
positions of the selected pair under investigation are in- 
dependently chosen at each check. In statistical physics, 
this situation resembles that of annealed disorder where 
random variables are allowed to evolve with time. 

The opposite of annealed disorder is quenched disor- 
der. Here, random variables do not evolve with time, 
i.e. they are quenched or frozen. This can be under- 
stood best for hard connectivity functions H{r) i.e. when 
77 = 00. For a particular realization of the Nj^ node posi- 
tions in sub-domain A^ the union of all connection regions 
IJj B{a.i) forms a cumulative connection region and grows 
with Nj( (see Fig. 3). Since nodes in B are uniformly dis- 
tributed, the probability of no bridging links is equal to 
(1 ~ I Ui B{s^^)\/V1s)^'^ , and so on average 



X = 1- 1- 



(IUg(a.)|).A 
Vb 



Nb 



(15) 



where (|lJj;B(ai)|)^ is the average area of the cumula- 
tive connection region due to all possible spatial config- 
urations of nodes in A. 

For small openings relative to the typical node connec- 
tion range tq, the cumulative connection region IJ^ B^&i) 
is a complicated object (see Fig. 3.b.) leaving little hope 
in obtaining an expression for its expected area in closed 
analytic form. Nevertheless, Eq. (15) provides further 



insight into the problem as it suggests that only nodes in 
one of the two sub-domains are quenched; the system is 
se?7ii-quenched. Therefore we may simplify (6) to obtain 
our main result 



Na Nb 

i=l 3 = 1 

,Nb, \^-^ 



l-{{{l-{X^3H.3)Bf')Ay 



(16) 
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where we recognize the term in the final brackets as the 
probability of node i not connecting with any of the N^ 
nodes in B. Note that Eq. (16) is symmetric under A <-> 
B and is also valid for higher dimensional domains. 



VII. COMPUTER SIMULATIONS AND 
DISCUSSION 

In this section we verify Eq. (16) through direct numer- 
ical integration using standard routines and comparing 
with computer simulations. The computer code written 
in C-l— I- uses Monte Carlo methods to estimate X for dif- 
ferent densities. At each run, an array of Ns randomly 
chosen coordinates in B is initiated. A random point is 
then generated in A, and is "checked" against all points 
in the array. The check is successful if the point in A has 
a line of sight and a connection link according to H{r) 
with at least one of the points in B. This is repeated 
up to iV_4 times, or until a successful check. Thus, the 
worst case computational cost is 0{N_^Nb). Note that an 
initial sorting of the array in ascending node-to-opening 
distance order, would typically offer a significant overall 
speed up due to the increased likelihood of a successful 
check. 

The simulation results are shown in Fig. 4 using solid 
markers for different values of ui = {0.01,0.05,0.1,0.2} 
using ?7 = 2, /3 = 1 and p_4 = pe in a square domain 



^ 



1.0 


- - ^^ £t t*j 


1 • \'_i^ it ■ Mt ■ ■ ■ 


0.8 
0.6 




■^ 


0.4 


'/'/y J^ 


w = 0.01 • 




nil 1 7 t^ 


u' = 0.05 - 


0.2 




w = 0.1 ♦ 


on 




IV = 0.2 - 








5 


10 


15 


1 


•M^ 


^yjj, , i^^^^^ 










^^^v*^'^-*^ ** 


^**i. 


*'*^**~*«~* 


0.1 




^s*V^''*^i^v, 




- ~?? 


'^. 




\ %\ 


V, 


t^^^ 


1 0.01 




vs 


V 


^^f_ 


! 1 




\ A^ 


> 


\, ^^ 


0.001 






N 




10-" 




\ 







p 



10 



FIG. 4: Comparison of numerical integration of (16) (solid 
curves), with computer simulations (solid markers), and the 
crude approximation of (12) (dashed curves) for different val- 
ues of w = {0.01,0.05,0.1,0.2} using t? = 2, /3 = 1 and 
Pa = Pb in a square domain as in Fig. 1. 



of side L = 10 as in Fig. 1. Numerical integration of 
Eq. (16) is shown using solid curves, and is seen to be 
in very good agreement with the computer simulations 
thus verifying the semi-quenched nature of the keyhole 
system set-up. 

The leading order prediction of Eq. (12) is shown in 
Fig. 4 using dashed curves. As expected, this approxi- 
mation is adequate at low densities p and small openings 
w, but worsens when increasing either p or w. Indeed, 
if the distance d.„ of the nearest node to the opening is 
much larger than w then d„ ~ p~^/^ and there is little to 
no overlap between connection regions (see Fig. 3.c.) and 
the independence assumption in Eq. (8) is reasonable. To 
estimate the critical density Pc at which (12) begins to 
diverge we set Pc = PA = Pb, and require 1 — X = 0(1), 
(e.g. X = e"i) giving pc ^ ^/W^w. 



VIII. CONCLUSION 

In this paper we have investigated the connectiv- 
ity properties of random networks in non-convex do- 
mains, and in particular when two networks in convex 
sub-domains may communicate through a small open- 
ing. We have shown that the connectivity properties 
of the complete network is a product of three indepen- 
dent events, the connectivity of each sub-network and 



X, the probability that a communication link exists be- 
tween them. Furthermore, we have shown that whilst 
both sub-networks have nodes randomly distributed in 
space, connectivity across the complete network exhibits 
semi-quenched disorder as only nodes in one of the two 
sub-domains are frozen random variables during the link 
formation process. We have benefited through this re- 
duction substantially, and have obtained general expres- 
sions for X (see Eq. (16)) and confirmed them through 
computer simulations. 

Further to the semi-quenched nature of the system, 
it seems clear that the keyhole problem is a challenging 
one due to the number of length scales involved. We have 
the typical size of system which is about ^V , the typical 
connectivity range tq ~ /3~^''', the typical distance be- 
tween neighbouring nodes in two dimensions p~^'^, and 
finally the width of the opening w. Hence, asymptotic 
methods requiring significant separation between each of 
the scales makes further analysis a formidable challenge. 
For example, even when w^tq there are a •priori three 
possible regimes: a) p"^/^<u;, b) w</9^^/^<ro, and c) 
p^^l"^ > rg. Note that in physical systems, more length 
scales may become relevant, e.g. the signal wavelength. 
Nevertheless, we have obtained leading order analytic 
expressions for X (see Eq. (12)) valid up to densities 
p,^./WJw. 

Familiar physically relevant settings where the above 
theory may be applicable are wireless networks in resi- 
dential, business or industrial buildings. Here, scattered 
devices form a machine to machine communication net- 
work whose functionality can be disrupted by the topo- 
logical features of the confining space, e.g. doorways 
and windows. Perhaps a more suitable application de- 
rives from ad hoc networks with minimal configuration 
requirements and quick deployment rates, often used for 
emergency situations like natural disasters or military 
conflicts. Here, the connectivity through small open- 
ings may correspond to that through narrow pathways, 
gaps and cracks encountered, e.g. in post earthquake 
urban environments, underground tunnels, or even in- 
side the human body [18]. Indeed, modern swarms of 
small robots can communicate, sense, collaborate, fetch 
and carry [19]. These technological advancements lend 
themselves to search and rescue operations in hazardous 
environments [20] and also to non-invasive medical diag- 
nostics and treatments by so called "nanobots" [21]. On 
the theoretical side, we expect our results to lead to fur- 
ther insight and understanding on network connectivity 
in non-convex domains. 
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